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Abstract We present deep HST/NICMOS observations peering through the outflow cavity of the protostellar candidate 
IRAS 04381+2540 in the Taurus Molecular Cloud-1. A young stellar object as central source, a jet and a very faint and close 
(0.6") companion are identified. The primary and the companion have similar colours, consistent with strong reddening. We 
argue that the companion is neither a shock-excited knot nor a background star. The colour/magnitude information predicts a 
substellar upper mass limit for the companion, but the final confirmation will require spectroscopic information. Because of its 
geometry, young age and its rare low-mass companion, this system is likely to provide a unique insight into the formation of 
brown dwarfs. 
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1. Introduction 
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The debat e on the origin of brown dwarfs has preceded their 
discovery (Reb olo etalJll995h : their small masses poses se- 
rious difficulty to a stellar-like formation scenario through 
the collapse of a molecular cloud core. Now, with hun- 
dreds of known brown dwarfs, the debate is still not set- 
tled: the most widely discussed formation scenarios in- 
clude grayJtati£nal_TOnapseJn highly turbulent cloud cores 
(e.g. IPadoan & Nordlun d||2004|) , early ejection from unstable 
multiple systems (Sterzik & Durisen 1998; Reirjurth & Clarke 
2001; Umbreit et al.|| 2004l) and photoevaporation of protostars 
lwMtworth^Zmnecke?l2004l) . Only few traces of the for- 
mation process survive the early evolution of these objects. 
Current studies focused on disk indicators, signs of accretion, 
binarity and velocity dispersion. We learned that brown dwarfs 
are often surrounded by warm circumstellar mat erial (e.g. 

I Javawardhana et all2003l) . with masses of ~Mj up jKlein et alJ 
2003) which they maintain over timescales similar to stars 
( Sterzik et al. 2004), and in one case it could be shown that the 



dust forms a disk (Pascucci et al. 2003). Because gravitational 
ejecti on leaves sufficient material in the inner disk dBate et alJ 
2002) to produce infrared excess emission and accretion signs, 
the ejection embryo hypothesis remains a possible formation 
scenario. The masses derive d from millimetre-wavelength ob- 
servations ( Kl ein et alJ20 03 ) are too uncertain to constrain pos- 
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sible disk truncation and therefore cannot exclude the ejection 
scenario. 

Direct observations of the youngest phases of substellar ob- 
jects, termed as proto-brown dwarfs, will be cruci al in evaluat- 
ing th e for mation scenarios. Re cent imaging by lHaisch et alJ 
(2004) and Duchene et al ] (12004 identified faint, possibly sub- 
stellar companion candidates to protostellar systems, but some 
of them might be background stars. To ensure the very young 
age of the proto-brown dwarfs, they should be observed within 
the massive protostellar envelope, possible only in exceptional 
cases. 

We present here observations of the first such sys- 
tem, IRAS 04381+2540, located in TMC-1 in the Taurus 
star-forming region with the discovery of a proto-brown 
dwarf candidate. Based on its spectral energy distribution 
and bolometric temperature, IRAS 04381+2540 has bee n 
classified as a Class I object (iHogerheiide & Sandelll EoOO). 
Ho wever, estimates of t he stellar mass (M = 0.3 ± 0.1M Q , 
IBrown & Chandlell Il999l) show tha t it is comparable to the 
envelope mass (M 



;t al.1 120031 indi- 



O.18M , 

eating that the objec t is still in its main accretion phase 
jMotte & Andrell200ll). IRAS 043 81+2540 drives a molecular 
outflow ( IBrown & Chandlerl !l999). with an inclinati on of i m 
50° (t he northern part tilted toward the observer. lchandler et alJ 
1996). Through the evacuated outflow cone we see into the 
dense molecular cloud core in which IRAS 0438 1+2540 is em- 
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Figure 1. High-resolution near-infrared view of 
IRAS 04381+2540 by HST/NICMOS in the broadband 
F160W filter. This image reveals a binary young stellar object 
(NIR A+B) and a collimated jet. 

bedded, making it an exceptional target for studies of very early 
stellar evolution. 

2. Observations and data reduction 

We used archival HST/NICMOS observations to study 
IRAS 04381+2540 at small spatial scales. The main parame- 
ters of the observations are summarised in Tab.[Q 

The data have been successfully pipeline-processed and re- 
quired only minor additional reduction work, which was car- 
ried out by simple IDL scripts. The images have been bad 
pixel filtered and the central columns bias-corrected in the stan- 
dard fashion. The final mosaic has been composed by cross- 
correlating the individual pointings. In order to enhance the 
fainter details around the central point source we subtracted 
synthetic point spread functions (P SF) calculated by using the 
Tiny Tim 6.1 jRrist & Hookl[l997i) . The scaling factors have 
been obtained by comparing the integrated counts in 3-pixel- 
radii apertures placed on the simulated and observed peaks. 

Aperture photometry has been performed on all images in 
an identical way. We used the IDL adaptation of the DAOPHOT 
routine to integrate the counts in 6-pixel-radius apertures cen- 
tered on the central point source. In order to obtain reliable 
photometry also for the fainter, secondary point source we car- 
ried out PSF photometry by subtracting simulated PSFs. By 
minimizing the subtraction residuals we measured the point 
source's brightness without significant contribution from the 
surrounding nebulosity. In the case of the F212N filter, due to 
the faintness of the secondary point source compared to the 
nebula, we could only estimate upper limits from the aperture 
photometry. Here, we also applied a second, small 3x3 pixel 
aperture to suppress the contamination from the nebula and 
thus obtained a stricter upper limit. 

3. Results 

As shown in Fig.^the NICMOS image reveals a detailed view 
of the reflection nebulosity including two point sources (IRAS 



0438 1+2540 NIR A and IRAS 04381+2540 NIR B) and a long 
optical jet at the heart of IRAS 04381+2540. In the following, 
we will refer to the brighter, primary point source as NIR A and 
to the fainter, secondary one as NIR B. NIR A is located in the 
apparent center of symmetry of the reflection nebula with the 
position of RA: 04'' 41 m 12'.'65 DEC: +25° 46' 35'/93 (J2000). 
This position is accurate within the spacecraft's pointing ac- 
curacy, typically better than 075. All NICMOS images show 
NIR A as a point source at a resolution of 0715. The secondary 
source NIR B is located at 0759 east and 07 15 north of NIR A. 
The source is unresolved in all three NICMOS images, while in 
the F160W image a possible marginal, faint extension is seen at 
the brightness level of the surrounding nebula, possibly a faint 
scattered-light envelope of NIR B. The photometry results are 
summarized in Tabled The deepest F160W image also reveals 
the presence of a ~3'/3 long, narrow line of intense emission, 
which we identified as a highly collimated jet. While the north- 
ern end of the jet vanishes in the background noise about 5" 
from NIR A, its southern end is overshined by the reflection 
nebula at about 1.2" from NIR A. The jet's direction points to 
NIR A identifying its driving source. 

4. Discussion 

The detection of the collimated jet proves that NIR A, its driv- 
ing source, is still in the early, accreting phase of its evolu- 
tion. By assuming an inclination of i = 50° for the jet its 
visible length is ~ 1000 AU. As often observed, the opti- 
cal jet in IRAS 04381+2540 is accompanied by a molecular 
outflow (Chandler et al. 1996). The present and former studies 
suggest the following picture: the young stellar object NIR A 
- at a transition stage in the Class 0/Class I boundary - is lo- 
cated close to the side of the dense molecular cloud core fac- 
ing the observer. The low-mass source NIR A drives a colli- 
mated jet and harbors a circumstellar accretion disk perpendic- 
ular to the jet's direction. The jet-powered outflow clears out 
a low-density cavity, whose wall becomes bright by scatter- 
ing the light of NIR A. The southern jet is inclined away from 
the observer and, by penetrating into the denser regions of the 
cloud, becomes strongly extincted and apparently fainter. The 
northern beam of the jet is inclined toward the observer and 
is visible inside the outflow cavity. The age of the system can 
be estimated on the basis of the ~ 10 5 yr typical lifetime of 
protostellar objects, consistent e.g. with th e chemical evolu- 
tion studies of the deuterium fractionation dSaito et aklkool) 
in IRAS 04381+2540. 

4. 1. The nature of the point sources 

The near-infrared colour information allows preliminary con- 
clusions on the nature of NIR A. When plotted in a colour- 
magnitude diagram (see Fig. the photometry of NIR A is 
consistent with that of a highly reddened (Ay « 43 mag) 
young stellar object. The mass estimate of ~0.6 M Q based on 
the compariso n of the near- infrared photometry to the 1 Myr 
isochrones of Bar affe et alJ dl998l) and is consistent with th e 
observed bolometric luminosity (0.7 L Q , IChandleretai1 l998i. 
A similar mass estimate was derived from gas kinematics 
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Instrument 


Filter 


Date 


Field of 


Exp. Time 


Flux of NIR A 


Flux of NIR B 


Note 






DD/MM/YY 


View 


[s] 


[/Jy] 


[mag] 


IpJy] 


[mag] 




HST/NICMOS 


F160W 


20/12/97 


20'xl9' 


1280 


511 


15.8 


38 


18.6 


1.4-1. 8/im 


HST/NICMOS 


F205W 


07/01/98 


19'xl9' 


256 


6295 


12.6 


363 


15.7 


1.75-2.35/im 


HST/NICMOS 


F212N 


07/01/98 


19'xl9' 


160 


7711 


12.3 


<1070 


>14.5 


2.121/iin, H2 


HST/NICMOS 


F212N 


07/01/98 


19'xl9' 


160 






<351 




0.13" aper. 



Table 1. Log of the infrared observations and results of the NICMOS photometry. The estimated photometric error is less than 
10% for all bands. In the F212N band only an upper limit can be reliably established due to the low contrast between the object 
and the background nebulosity. Here, the use of a smaller aperture reduces the contamination from the nebulosity and places a 
stricter upper limit. The F160W filter also includes [Fe II] line emission at 1.644 pm. 



04381 +2540 NIR A 




1 60W-F205W [n- 



Figure 2. Colour-magnitude diagram based on HST/NICMOS 
photometry for the IRAS 043 81+2540 NIR A and B sources. 
The 1 Myr-old isochrone from Baraffe et alJ (119981) is overplot- 
ted. The numbers indicate the masses in solar units. The inter- 
stellar reddening vectors for Ay=43 mag and 37 mag are also 
plotted. Thick lines mark the approximate mass range as con- 
strained by the photometric errors. 



(0.3±0.1 Mq, iBrown & Chandlerlll999l) . The case of NIR B 
is more compelling; as NIR B does not have any evident link 
to its environment suggesting it to be a companion of NIR A, 
it could also be an unrelated background/foreground star, or 
an excited knot in the cavity wall. These three possibilities are 
briefly discussed below. 

(1) Foreground/background source: The high reddening of 
NIR B (Ay ~37 mag) excludes the possibil ity of an un re- 
lated foreground star. The comparison (B alazs et al.ll2004T) of 
the lWainscoatet all \ 19921) model to the 2MASS Point Source 
Catalogue predicts A\j «16 mag extinction towards our source 
and the probability of a chance background star appearing 
within ~1.5 square arcsecond to be « 10~ 4 . The additional 
A v >20 mag extinction observed towards IRAS 04381+2540 
makes such a coincidence even less likely. 

(2) A shock-excited knot associated to the jet: Although 
several such examp les are known, e.g., in L1551 NE (e.g. 
iReiourth et aljEoOOl) . the source NIR B differs from these ob- 
jects in three important aspects: it is visible and unresolved in 
the broad-band images; and its brightness in the F212N filter 
is dominated by continuum emission (i.e. no shocked H2 emis- 
sion is visible). 

(3) Companion to NIR A: NIR B has near-infrared colours very 
similar to NIR A; a compariso n to the 1 Myr-old substellar 
models of Baraffe et suggests a reddening almost 
identical to NIR A (see Fig. 0. In addition to the matching 



colours, the very proximity of NIR A and B further reinforces 
their physical link. The alternative scenarios being very im- 
probable, in the following we tentatively accept NIR B to be 
a point source related to NIR A. 

4.2. Probable substellar mass for NIR B 

The comparison of NIR B brightness and colour (see Fig. 
argues for a substellar mass of Mb ~ 50 ± 4 Mj up . The error 
bar represents the photometric error, but in the following we 
discuss several other factors which influence this estimate. 

(1 ) Isochrones: Our object is likely to be younger than the 
youngest available isochrones dating 1 Myr. The yet un- 
known initial conditions of brown dwarf formation probably 
introduce sign ificant uncertainities even for the 1 Myr tracks 
(Bar affe et alJl2002l) . Comparing the luminosity and effective 
temper ature in the 40 -60 Mj up range to the 1 Myr-old mod- 
els by iBurrows et alJ (I19971) . we find that the latter models 
associate 3-5 Mj up smaller masses to the given luminosi- 
ties, suggesting an even smaller mass for NIR B. Although 
no isochrones of the appropriate age are yet available, using 
overaged isochrones lead to the overestimation of the mass and 
therefore enables us to set an upper mass limit for NIR B of 
-50 M Jup . 

(2) Reddening: The reddening law plays a crucial role when 
interpreting objects as heavil y reddened a s NIR B. In Fig.|2]we 
used the reddening law by Mathis ( 1990) and now we explore 
how the estimate varies with the reddening law's slope. A study 
of the reddening toward the Taurus region by Whittet et alJ 
(2001) showed that in the direction of dense clouds Ry= 3.5 
- 4.0 and might reach 4.5 for the densest regions of the cloud. 
Disentangling the contribution to NIR B's reddening by the dif- 
ferent dust types is not possible. By making the worst-case es- 
timate (assuming that all the extinction is caused by dust with 
Ry as high as 4.5) would modify NIR B's upper mass limit to 
~90 Mj up , i.e. slightly above the substellar limit. We conclude, 
that although the reddening law influences the mass estimates, 
it is very unlikely to shift the upper mass limit into the stellar 
regime. 

(3) Anisotropic emission: The emission of embedded stars is 
often influenced by the circumstellar environment via obscu- 
ration and lig ht scatteri ng, exemplified by edge-on disks (see, 
e.g. lPadgett et alJl9 99). In these cases both the colour and the 
magnitude of the objects are altering mass estimates on the ba- 
sis of theoretical isochrones. Circumstellar material obscuring 
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more than 66% of NIR B's emission could apparently shift the 
luminosity of a very low mass star to substellar luminosities. 
Alternatively, NIR B's emission might originate partly by its 
envelope, scattering the light of NIR A; then NIR B's mass 
might be even smaller than estimated. Entirely excluding these 
possibilities is impossible without detailed spectroscopic and 
polarimetric follow-up observations but we argue that such a 
strong scattered-light nebula would li kely be r esolvable by our 
observations. While a study by Padg ett et alJ 0999) shows 7 
such nebulae around edge-on disks to have typical sizes of 
~800 AU, NIR B is an unresolved point source at our ~20 AU 
resolution. Thus, should NIR B be a reflection nebula, its point 
source nature would be rather exceptional. 

4.3. NIR B as a brown dwarf in formation 

The concept of NIR B being a substellar companion of NIR A 
opens up exciting possibilities. Although several stars with 
brown dwarf companions are known, this object is about an or- 
der of magnitude younger than the previously known youngest 
brown dwarf com panion to the T Tauri star GG Tau Bb 
ftVhite et alj F999). The presence of the optical jet confirms 
active accretion in the system and suggests that NIR B is a 
brown dwarf in formation. An important question is whether 
NIR B will be able to accrete sufficient mass to become a star 
during the future accretion processes. Observations of the B1-7 
line can determine the accretion rates of NIR A and NIR B, 
and show whether the current mass difference will increase 
(to produce a low-mass star with brown dwarf companion) or 
decrease (toward an equal-mass binary situation). Additional, 
longer wavelength high-resolution observations can identify 
whether NIR B is embedded in the disk of NIR A or harbors 
its own accretion disk, an important clue on the formation his- 
tory of the system. Should NIR B be in the disk, it is likely to 
gravitationally interact with the disk and its mass is probably 
sufficient to induce observable structures, similar to the ones 
predicted for massive planets embedded in Solar- ty pe stars' 
disks (e.g. lBrvden et al)l200Cl:lb'Angelo et alJl2003l) allowing 
the observational test of these models. 

5. Conclusions 

Based on our HST/NICMOS infrared imaging the main con- 
clusions of this work are the following: 

(i) The central source of the IRAS 04381+2540 system is iden- 
tified as a young, low-mass binary stellar object. 

(ii) The accreting primary component drives a collimated jet. 

(iii) The secondary component is probably the youngest known 
brown dwarf so far, providing a unique opportunity to study the 
formation of substellar objects. 
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